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Abstract 

This  paper  is  on  beamforming  in  wireless  relay  networks  with  perfect  channel  information 
at  relays,  the  receiver,  and  the  transmitter  if  there  is  a  direct  link  between  the  transmitter  and 
receiver.  It  is  assumed  that  every  node  in  the  network  has  its  own  power  constraint.  A  two-step 
amplify-and-forward  protocol  is  used,  in  which  the  transmitter  and  relays  not  only  use  match 
filters  to  form  a  beam  at  the  receiver  but  also  adaptively  adjust  their  transmit  powers  according 
to  the  channel  strength  information.  For  a  network  with  any  number  of  relays  and  no  direct  link, 
the  optimal  power  control  is  solved  analytically.  The  complexity  of  finding  the  exact  solution 
is  linear  in  the  number  of  relays.  Our  results  show  that  the  transmitter  should  always  use  its 
maximal  power  and  the  optimal  power  used  at  a  relay  is  not  a  binary  function.  It  can  take  any 
value  between  zero  and  its  maximum  transmit  power.  Also,  surprisingly,  this  value  depends  on 
the  quality  of  all  other  channels  in  addition  to  the  relay’s  own  channels.  Despite  this  coupling 
fact,  distributive  strategies  are  proposed  in  which,  with  the  aid  of  a  low-rate  broadcast  from  the 
receiver,  a  relay  needs  only  its  own  channel  information  to  implement  the  optimal  power  control. 

*This  work  was  supported  in  part  by  ARO  under  the  Multi-University  Research  Initiative  (MURI)  grant  #W91 1NF- 
04-1-0224. 
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Simulated  performance  shows  that  network  beamforming  achieves  the  maximal  diversity  and 
outperforms  other  existing  schemes. 

Then,  beamforming  in  networks  with  a  direct  link  are  considered.  We  show  that  when  the 
direct  link  exists  during  the  first  step  only,  the  optimal  power  control  at  the  transmitter  and 
relays  is  the  same  as  that  of  networks  with  no  direct  link.  For  networks  with  a  direct  link  during 
the  second  step  only  and  both  steps,  recursive  numerical  algorithms  are  proposed  to  solve  the 
power  control  problem.  Simulation  shows  that  by  adjusting  the  transmitter  and  relays’  powers 
adaptively,  network  performance  is  significantly  improved. 


1  Introduction 

It  is  well-known  that  due  to  the  fading  effect,  the  transmission  over  wireless  channels  suffers  from 
severe  attenuation  in  signal  strength.  Performance  of  wireless  communication  is  much  worse  than 
that  of  wired  communication.  For  the  simplest  point-to-point  communication  system,  which  is  com¬ 
posed  of  one  transmitter  and  one  receiver  only,  the  use  of  multiple  antennas  can  improve  the  capacity 
and  reliability.  Space-time  coding  and  beamforming  are  among  the  most  successful  techniques  de¬ 
veloped  for  multiple- antenna  systems  during  the  last  decades  [1,2].  However,  in  many  situations, 
due  to  the  limited  size  and  processing  power,  it  is  not  practical  for  some  users,  especially  small 
wireless  mobile  devices,  to  implement  multiple  antennas.  Thus,  recently,  wireless  network  commu¬ 
nication  is  attracting  more  and  more  attention.  A  large  amount  of  effort  has  been  given  to  improve 
the  communication  by  having  different  users  in  a  network  cooperate.  This  improvement  is  conven¬ 
tionally  addressed  as  cooperative  diversity  and  the  techniques  cooperative  schemes. 

Many  cooperative  schemes  have  been  proposed  in  literature  [3-21].  Some  assume  channel  in¬ 
formation  at  the  receiver  but  not  the  transmitter  and  relays,  for  example,  the  noncoherent  amplify- 
and-forward  protocol  in  [8,  9]  and  distributed  space-time  coding  in  [10].  Some  assume  channel 
information  at  the  receiving  side  of  each  transmission,  for  example,  the  decode-and-forward  proto¬ 
col  in  [8, 12]  and  the  coded-cooperation  in  [13].  Some  assume  no  channel  information  at  any  node, 
for  example,  the  differential  transmission  schemes  proposed  independently  in  [14-16].  The  coher- 
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ent  amplify-and-forward  scheme  in  [9, 11]  assumes  full  channel  information  at  both  relays  and  the 
receiver.  But  only  channel  direction  information  is  used  at  relays.  In  all  these  cooperative  schemes, 
the  relays  always  cooperate  on  their  highest  powers.  None  of  the  above  pioneer  work  allow  relays 
to  adjust  their  transmit  powers  adaptively  according  to  channel  magnitude  information,  and  this  is 
exactly  the  concern  of  this  paper. 

There  have  been  several  papers  on  relay  networks  with  adaptive  power  control.  In  [22,  23], 
outage  capacity  of  networks  with  a  single  relay  and  perfect  channel  information  at  all  nodes  were 
analyzed.  Both  work  assume  a  total  power  constraint  on  the  relay  and  the  transmitter.  A  decode-and- 
forward  protocol  is  used  at  the  relay,  which  results  in  a  binary  power  allocation  between  the  relay 
and  the  transmitter.  In  [24],  performance  of  networks  with  multiple  amplify-and-forward  relays  and 
an  aggregate  power  constraint  was  analyzed.  A  distributive  scheme  for  the  optimal  power  allocation 
is  proposed,  in  which  each  relay  only  needs  to  know  its  own  channels  and  a  real  number  that  can 
be  broadcasted  by  the  receiver.  Another  related  work  on  networks  with  one  and  two  amplify-and- 
forward  relays  can  be  found  in  [25].  In  [26],  outage  minimization  of  single-relay  networks  with 
limited  channel-information  feedback  is  performed.  It  is  assumed  that  there  is  a  long-term  power 
constraint  on  the  total  power  of  the  transmitter  and  the  relay.  In  this  paper,  we  consider  networks 
with  a  general  number  of  amplify-and-forward  relays  and  we  assume  a  separate  power  constraint  on 
each  relay  and  the  transmitter.  Due  to  the  difference  in  the  power  assumptions,  compared  to  [24], 
analysis  of  this  new  model  is  more  difficult  and  totally  different  results  are  obtained. 

For  multiple-antenna  systems,  when  there  is  no  channel  information  at  the  transmitter,  space- 
time  coding  can  achieve  full  diversity  [1].  If  the  transmitter  has  perfect  or  partial  channel  infor¬ 
mation,  performance  can  be  further  improved  through  beamforming  since  it  takes  advantage  of 
the  channel  information  (both  direction  and  strength)  at  the  transmit  side  to  obtain  higher  receive 
SNR  [2] .  With  perfect  channel  information  or  high  quality  channel  information  feedback  from  the 
receiver  at  the  transmitter,  one-dimensional  beamforming  is  proved  optimal  [2,27,28].  The  more 
practical  multiple- antenna  systems  with  partial  channel  information  at  the  transmitter,  channel  statis¬ 
tics  or  quantized  instantaneous  channel  information,  are  also  analyzed  extensively  [29-33].  In  many 
situations,  appropriate  combination  of  beamforming  and  space-time  coding  outperforms  either  one 


3 


of  the  two  schemes  alone  [34-37].  In  this  paper,  we  will  see  similar  performance  improvement  in 
networks  using  network  beamforming  over  distributed  space-time  coding  and  other  existing  schemes 
such  as  best-relay  selection  and  coherent  amplify-and-forward. 

We  consider  networks  with  one  pair  of  transmitter  and  receiver  but  multiple  relays.  The  receiver 
knows  all  channels  and  every  relay  knows  its  own  channels  perfectly.  In  networks  with  a  direct 
link  (DL)  between  the  transmitter  and  the  receiver,  we  also  assume  that  the  transmitter  knows  the 
DL  fully.  A  two-step  amplify-and-forward  protocol  is  used,  where  in  the  first  step,  the  transmitter 
sends  information  and  in  the  second  step,  the  transmitter  and  relays,  if  there  is  a  DL,  transmit.  We 
first  solve  the  power  control  problem  for  networks  with  no  DL  analytically.  The  exact  solution  can 
be  obtained  with  a  complexity  that  is  linear  in  the  number  of  relays.  Then,  to  perform  network 
beamforming,  we  propose  two  distributive  strategies  in  which  a  relay  needs  only  its  own  channel 
information  and  a  low-rate  broadcast  from  the  receiver.  Simulation  shows  that  the  optimal  power 
control  or  network  beamforming  outperforms  other  existing  schemes.  We  then  consider  networks 
with  a  DL  during  the  first  transmission  step,  the  second  transmission  step,  and  both.  For  the  first 
case,  the  power  control  problem  is  proved  to  be  the  same  as  the  one  in  networks  without  the  DL.  For 
the  other  two  cases,  recursive  numerical  algorithms  are  provided.  Simulation  shows  that  they  have 
much  better  performance  compared  to  networks  without  power  control.  We  should  clarify  that  only 
amplify-and-forward  is  considered  here.  For  decode-and-forward,  the  result  may  be  different  and  it 
depends  on  the  details  of  the  coding  schemes. 

The  paper  is  organized  as  follows.  In  the  next  section,  the  relay  network  model  and  the  main 
problem  are  introduced.  Section  3  works  on  the  power  control  problem  in  relay  networks  with  no  DL 
and  Section  4  considers  networks  with  a  DL.  Section  5  contains  the  conclusion  and  several  future 
directions. 


2  Wireless  Relay  Network  Model  and  Problem  Statement 

Consider  a  relay  network  with  one  transmit-and-receive  pair  and  R  relays  as  depicted  in  Fig.  1. 
Every  relay  has  only  one  single  antenna  which  can  be  used  for  both  transmission  and  reception. 
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Denote  the  channel  from  the  transmitter  to  the  ith  relay  as  /,  and  the  channel  from  the  zth  relay 
to  the  receiver  as  gt.  If  the  DL  between  the  transmitter  and  the  receiver  exists,  we  denote  it  as 
/0.  We  assume  that  the  transmitter  knows  /0,  the  zth  relay  knows  its  own  channels  /,;  and  g, ,  and 
the  receiver  knows  all  channels  /0 ,  ,/'i ,  •  •  • ,  frt  and  gi, . . . ,  qr.  The  channels  can  have  both  fading 
and  path-loss  effects.  Actually,  our  results  are  valid  for  any  channel  statistics.  We  assume  that  for 
each  transmittsion,  the  powers  used  at  the  transmitter  and  the  it h  relay  are  no  larger  than  P0  and  Pt, 
respectively.  Note  that  in  this  paper,  only  short-term  power  constraint  is  considered,  that  is,  there  is 
an  upper  bound  on  the  average  transmit  power  of  each  node  for  each  transmission.  A  node  cannot 
save  its  power  to  favor  transmissions  with  better  channel  realizations. 

We  use  a  two-step  amplify-and-forward  protocol.  During  the  first  step,  the  transmitter  sends 
a(y\fp)8.  The  information  symbol  s  is  selected  randomly  from  the  codebook  S.  If  we  normalize  it 
as  E|s|2  =  1,  the  average  power  used  at  the  transmitter  is  a(2 P0 .  The  zth  relay  and  the  receiver,  if  a 
DL  exists  during  the  first  step,  receive 

T'i  =  ao^PofiS  +  Vi  and  xx  =  ao^/Pofos  +  w1:  (1) 


respectively.  Vi  and  tig  are  the  noises  at  the  zth  relay  and  the  receiver  at  Step  1.  We  assume  that  they 
are  CAf(().  1).  During  the  second  step,  the  transmitter  sends  3\y\fp)Cj0(i s,  if  a  DL  exists  during  this 
step.  At  the  same  time,  the  zth  relay  sends 


ti  Oti 


P 


l  +  «g|/,|2P0 


The  average  transmit  power  of  the  zth  relay  can  be  calculated  to  be  a] Pi.  If  we  assume  that  /0  keeps 
constant  for  the  two  steps,  the  receiver  gets 

R 

X2=f3oy/Pofoejd°s  +  ^  gtU  +  w2 

i=  1 


R 


—  \/Pq  (  Pofo^e°  +  OL 0 


aifigi£j()l\[Pi 


\  ^  oiigie^1  \fPi 

>  — ,  Vj  +  W2-  (2) 

Vl  +  al\M2po 


v  tr  n/i+«oI  M2po/ 

w2  is  the  noise  at  the  receiver  at  Step  2,  which  is  also  assumed  to  be  CJf( 0, 1).  Note  that  if  the 
transmitter  sends  during  both  steps,  we  assume  that  the  total  average  power  it  uses  is  no  larger  than 
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P0.  With  this,  the  total  average  power  in  transmitting  one  symbol  is  no  larger  than  -  Clearly, 

the  coefficients  a0,  oti, ,  aR  are  introduced  in  the  model  for  power  control.  The  power  constraints 
at  the  transmitter  and  relays  require  that  a  1  +  <1  and  0  <  a,  <  1. 

Our  network  beamforming  design  is  thus  the  design  of  90,  6 1,  •  •  •  ,9R  and  a0)  ,  aR, 

such  that  the  error  rate  of  the  network  is  the  smallest.  This  is  equivalent  to  maximize  the  receive 
SNR,  or  the  total  receive  SNR  of  both  branches  if  a  DL  exists  during  the  first  step.  From  (2),  we 
can  easily  prove  that  an  optimal  choice  of  the  angles  are  6*0  =  —  arg  /0  and  0,  =  —  (arg  f  ,  +  arg  cji). 
That  is,  match  filters  should  be  used  at  relays  and  the  transmitter  during  the  second  step  to  cancel 
the  phases  of  their  channels  and  form  a  beam  at  the  receiver.  We  thus  have 


—  /~P  ( R  \  f  \  L  Y^  ai\fi9i\VPi 

%2  —  V  Pq  I  Po  l/o  |  +  «0  2_^ 


R 


,  , _ _ ,  s + v  -skMS=e-i^f.Vi + W2.  (3) 

V  i= 1  \/l  +  al\fi\2P0  J  i=1  y/l  +  al\fi\2P0 

What  is  left  is  the  optimal  power  control,  i.e.,  the  choice  of  a0>  Po,  a i,  •  •  • ,  otR.  This  is  also  the  main 
contribution  of  our  work. 


3  Optimal  Relay  Power  Control 

In  this  section,  we  investigate  the  optimal  adaptive  power  control  at  the  transmitter  and  relays  in  net¬ 
works  without  a  DL.  Section  3.1  presents  the  analytical  power  control  result.  Section  3.2  comments 
on  the  result  and  gives  distributive  schemes  for  the  optimal  power  control.  Section  3.3  provides 
simulated  performance. 


3.1  Analytical  Result 


With  no  DL,  we  have  /30  =  0  and  x \  =  0.  From  (3),  the  receive  SNR  can  be  calculated  to  be 

.  2 

'R  oii\figi\-/Pi 
v/l+ogl/iPPo 


«oPo  Z)i=i 


1  I  ai\9i\2Pi 


6 


It  is  an  increasing  function  of  «o.  Therefore,  the  transmitter  should  always  use  its  maximal  power, 
i.e.,  a*  =  1.  The  receive  SNR  is  thus: 


p  |  y^R  ai\ftgt\'/Fri  \ 

0  \^i=1  Vl+\fi\2Po  ) 

1  ,  y^R  a1  \m\2Pi 
1  "T  2^i= l  i+iyPPo 

Before  going  into  details  of  the  SNR  optimization,  we  first  introduce  some  notation  to  help 
the  presentation.  (•,  •)  indicates  the  inner  product.  ||  ■  ||  indicates  the  2-norm.  P  indicates  the 
probability,  at  denotes  the  /th  coordinate  of  vector  a  and  denotes  the  /c-dimcnsional  vector 

r  i  T 


a?|  •  •  •  aik  ,  where  -T  represents  the  transpose.  If  a,  b  are  two  A-dimcnsional  vectors, 

a  ^  b  means  a j  <  bi  for  all  i  =  1 , ,R.  Op  is  the  72-dimensional  vector  with  all  zero  entries. 
Denote  the  set  Op  ^  y  ^  a  or  equivalently,  0  <  yi  <  a,i  for  i  =  1, ...,/?,  as  A.  For  1  <  k  <  R  —  1, 
denote  the  set  0fc  r<  yh,...,ik  A  a as  A where  {A, ...  ,4}  is  a  fc-subset  of  {1, ... ,  72}. 


Define 


\fl9l\VP~i 

\9lWP~l 

X  = 

0^1 

i  b  = 

vA+IM2^ 

,a  = 

x/i+l/iPRo 

,  and  A  =  diag(a), 

°lr 

I/r9rI\/Rr 

\9rWPr 

L  V1+i/«i2po  -J 

L  v^+l/flPRo  J 

where  diag(a)  indicates  the  diagonal  matrix  whose  zth  diagonal  entry  is  a,.  With  the  transformation 
y  =  Ax,  or  equivalently,  x  =  A-1  y,  we  have 


SNR 


P  (b^)2 
°1  +  ||Ax||2 


Po 


(c,y)2 
1  +  llyll2’ 


where 


y/l+\fi\2Pb 

\9i\VK 

0 

\fl9l\VPl 

\/i+l/il2Ro 

"  lAI  " 

0 

a/i+I/rPRo 

\9rWPr 

\fR9R\VPR 

-  A/1+Pi?|2p’°  - 

/r 

The  receive  SNR  optimization  problem  is  thus  equivalent  to 


max  AWT  y€A. 

y  1  +  ||y||2 


(4) 
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The  difficulty  of  the  problem  lies  in  the  shape  of  the  feasible  set.  If  y  is  constrained  on  a  hypersphere, 


that  is,  ||y||  =  r,  the  solution  is  obvious  at  least  geometrically.  Given  that  ||y||  =  r, 


c,y  2  r2  c  2  2 

,  ,  II  ||2  =  cos  V’ 

1  +  ||y|r  1  +r2 

where  tp  is  the  angle  between  c  and  y.  The  optimal  solution  should  be  the  vector  which  has  the 
smallest  angle  with  c.  Thus,  we  decompose  (4)  as 

2 


max 


r  l+r2V||y||=r 


I  max  (c,  y)  )  s.t.  y  G  A  and  0  <  r  <  ||a| 


(5) 


Since  P(aj  >  0)  =  1  and  P(cj  >  0)  =  1,  we  assume  that  a*  >  0  and  ct  >  0.  Define 


•,  o,.,,  ;; 

a3  I 9j\VP3 


Order  o,  as 


(6) 


Vi  ^  0T2  ^  ‘ 


(7) 


(ri,  t2, . . . ,  tr)  is  thus  an  ordering  of  (1,2,...,  i?).  Define 


P)=0, 


r1=0T11||c||  =  Jd)-2\\c 


T\  II  T2viTR\ 


+  a2  , 


r2= \/d>^2\\c. 


II 2  _i_  q2  — — 

T2  11^7-2, . ..,TR  II  I 


\ 


^r22|lCT3v,'ri?  | 


i—  1 


rR-i= 


rR= 


\ 


lip  I 

Rr-1 


Ft- 2 

E 

j=i 


a2  = 


-R-l 


A  ^ 


-i  r-'Tfl  i 


2=1 


\ 


-R-l 


^rR\CrR\2  +  J2 


a2.  =  Hall. 


2=1 


Since  oT]  ,  >  cf)Tj,  we  have  fj_  |  <  r3  for  j  =  1, . . . ,  R.  Thus,  the  feasible  interval  of  the  radius, 
[0,  ||a||],  can  be  decomposed  into  the  following  R  intervals: 


[0,  || a||]  =  [r0,ri]  U  [ri,r2]  U  ■  ■  ■  U  [rR-2,rR-i\  U  [rR_i,rR]. 
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We  denote  T,:  =  [ri:  ri+ 1]  for  i  —  0, . . . ,  R  —  1.  Thus,  (5)  is  equivalent  to 

2 


max  max 


i=i,...,R  rer\  1  +  r2 


max  (c,  y) 
||=rer;,yeA 


We  have  decomposed  the  optimization  problem  into  R  subproblems.  We  now  work  on  the  zth 
subproblem: 

1  /  \ 2 

(8) 


max 


ref;  1  +  r2 

Denote  the  solution  of  the  inner  optimization  problem, 


max  (c,  y) 
=reri)yeA  '  7  ' 


max  (c,  y), 

|!=reri,yeA 


(9) 


as  zW  We  have  the  following  two  lemmas. 
Lemma  1.  z-  =  a,  for  j  =  T]. ... .  r,. 


Proof.  We  prove  this  lemma  by  contradiction.  Assume  that  zf  <  dj  for  some  j  G  {n, . . . ,  r,}.  We 


first  show  that  there  exists  an  l  G  {ri+i, . . . 

,  rR}  such  that  -f- 

m  G  {ri+i, . . . ,  rR}.  We  have 

\  zw 

\  c  — — 
_  r . 

cj 

<  Cmg  = 

i 

\  2  /  ...  x  2 

lIzWll  - 
l|z  11  “  \ 

Y 

j  +  [Zt™) 

\ 

m=  1 

m=i-\- 1 

i 

R 

<  \ 

Y  <,  +  S  0/ 

\ 

ra=l 

m=i+ 1 

=  \ 

fj  |  CTi+l 

,..^ii2  +  E< 

ra=l 

^  \ 

0Ti2|  Cri+i 

ra=l 

(i) 


becasue  of  (7) 


=  r. 


Z<P  ZW 

This  contradicts  ||zw||  g  Tj.  thus,  there  exists  an  /  G  {ri+ 1, . . . ,  rR}  such  that  —  < 


Define  another  vector  z'  as  z'„  =  +  5,  z,  = 


,(*) 


28 z^  -  82,  and  z'  =  Zm  for 


W 


where 


c2\  -1  /„«  ZW' 


0  <  <5  <  min  <(  2cj  (  1  +  -^  )  (  -  f-  )  , 


(i)  (i) 

~  Zji  aj  ~  zj 
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Since  we  have  assumed  that  >?  <  a3  and  have  just  proved  that  such  5  is  achievable.  To 

contradict  the  assumption  that  is  the  optimal,  it  is  enough  to  prove  the  following  two  items: 


zV  JO 


1.  z'  is  a  feasible  point:  ||z'||  =  r  and  z'  e  A, 

2.  (c,z«)  <  (c,z'). 


From  the  definition  of  z ',  we  have 


|z'||2  =  (z'V+u; 


/)2+  ~  +  (*7 

m^j,l 


(0 


-2fcf-52+E 


=  ||z»||2  =  r2. 


m/j,Z 


Since  0  <  5  <  a,-  —  we  have  0  <z'  <  cij.  Also,  since  0<c5<cW(zj',j  T  ^  z^ 


JO 


(0 


(i) 

z)  ,  we 


can  easily  prove  that  z[  = 


JO 


—  25zf>  —  52  >  0  and  z\  <  z)v  <  ai.  Thus,  z'  e  A.  The  first 


JO 


JO 


c2\-l  /2W 

item  has  been  proved.  For  the  second  item,  since  5  <2Cj[l  +  ^J  I  -L-  -  ] ,  we  have 


1  +  %)  S2  <  2c,  I  ^  )  5  =  2  {  -  zf  ]  5 


JO 


■'i 

2  C 


+  3*2-2^<  K 


Q 


AO; 


Q 


Q 


JO 


2zf  5  -  52  =  zf 


QZj  -CjS  <  ciz[ 

Ciz[l)  +  c,-zjl)  -  (cjZj  +  Cj-z')  <  0 
(c,zw)  <  (c,z). 


□ 


Lemma  2.  z} ;  = 


(0  _  \/^2-ELi< 


Ti+1>  — )TflH 


■Cj/or  j  =  Ti+ 1, . . .  ,tr. 


(i) 

Proof.  From  Lemma  1,  z  =  a3  for  j  -=  T\. ... .  rt.  Thus,  (9)  can  be  written  as 


,,  .  E  aTmCTm  +  (C 

li=rer;,yeA 

m=l 


yri+i,...,rfl) 


E6t-  + 


max 


m=l 


l=Vr2-^m=l  “?m> 


yrl+i,...  }  • 


i’»Tj+1,...,T)j6ATj+lv..,-rfl 
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Define  A  = 


fTi+i,...,TR  | 


A  =  !t  is  obvious  that  (cri+1,...,rii,yri+1,...,Tfl)  <  (cTi+lr..,rjl,  AcTi+li...iTfl)  for  all 


*+l 


=  \  r 


-  .  U-2 


XTm=i  arm-  *n  other  words,  to  maximize  the  inner  product,  yTi+lr..>Tfl  should 


have  the  same  direction  as  cT.+1  rjr  Thus,  we  only  need  to  show  that  this  direction  is  feasible  for 
r  e  Tj.  This  is  equivalent  to  show  that  AcTi+1)...iTB  6  ATi+1).„iTH  for  any  r  e  T*.  We  can  easily  prove 
that 

r  e  Ti  A  6  Oj, 


where  Q,  = 


for  i  =  0, . . . ,  R  —  1.  Thus,  for  any  r  e  T,  and  j  =  ri+ 1, . . . ,  rR,  we  have 


0  <  \cj  <  0r.+iCj  <  (j)j  Cj  cij.  Hence,  AcTi+1)...jTH  €  AT4+1j...jTh. 


□ 


Combining  Lemma  1  and  Lemma  2,  we  have 


(i) 

= 


Ac,- 


J=TU...,Ti 
j  ^i+1?  ■  •  •  j  7 R 


(10) 


and  thus 


max  (c,  y)  =  V  bTm  +  A||cTi+li...,TB| 

||=reri,yeA  “ 

m=  1 


(ll) 


We  have  solved  the  inner  optimization  of  Subproblem  i.  The  solution  of  the  R  subproblems  can 
thus  be  obtained. 


Lemma  3.  The  solution  of  Subproblem  0  is  y(0)  =  0  1  c.  For  1  <  i  <  R  —  1,  the  solution  of 
Subproblem  i  is  yit}  that  is  defined  as 


mm 


r  i+su 

l  eUh 


-1  “t 
L  rm 


6"1 

r<+i 


j  =  n, . . . ,  n 
j  •  •  •  5  Tft* 


(12) 


Proof.  From  (11),  Subproblem  i  is  equivalent  to  the  following  1-dimensional  optimization  problem: 


max 


Era=l  +  ||cri+i,...,rH| 


2A 


AeOi  1  +  Y!m=i  +  ||cTi+i,...,-rJ|2A2 


(13) 


When  i  =  0,  (13)  is  equivalent  to  max^n* 


||c||4A2 


*  l+||c||2A2 


.  Since 


||c||4A2 
1-|- 1|  c|  1 2  A2 


is  an  increasing  function  of 


A,  its  maximum  is  at  A  = 


L  —  l 
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For  i  —  1, . . . ,  R  —  1,  Define 


6(A)  = 


y*  1  bT  +  lie 

Z— /m=l  Tm  1  II 


Ti+l,...,TR  | 


2A 


1  +  ELl  arm  +  llCri+i,...,rH||2A2 


We  have, 


^  (Em=1  bTm  +  ||cTi+1)...jTH||  AJ  ||cTi+lj... 

fi  +  EL.  <  +  lk,+> . jw 


TR  I 


1  +  ^  “L  ~~  6mA 


771—1 


771=1 


^  D  if  \  ^  1  +  EA  =  1  “r,„  Qrlr|  9?  ^  n  if  \  ^  !  +  Em  =  l  aT,„  .  gQ,  j+E m  =  l  arm  <  ^,-1 


Thus,  ^  >  0  if  A  < 


and  ^  <  0  if  A  > 

Em=l  dA  Em=l  bTm 


EUl  b-r 


—  rri+i’ 


1+V*  a2 

the  optimal  solution  is  reached  at  A  =  —  m=1  Tm 


Em=l 

A  =  0“1i.  From  (10),  Subproblem  i  is  solved  at  y6)  as  defined  in  (12). 


.  Otherwise,  the  optimal  solution  is  reached  at 

□ 


Now,  we  can  work  on  the  relay  power  control  problem  presented  in  (4). 


(b,x«)2 


Theorem  1.  The  solution  of  the  SNR  optimization  is 

x*  =  arg  max  - 

^i=i  R-i,  1  +  IlylxWH2’ 

where  for  i  —  1, . . . ,  R  —  2,  x^)  is  defined  as 


(14) 


(i) 

Xj  = 


i+e; 


V*  h 

2^m= 1  ®Tm 


J  =  Ti,...  ,Ti 
j  Tj+1  :•••■>  Tfi 


and  the  binary  function  I,  is  defined  as 


h  = 


0  if 


if  1  +  Em=l  “2 


E  7n=l  ^Tm 

1  otherwise 


—  >  0”1 


(15) 


(16) 


JR-l) 


is  defined  as 


x 


(R-l) 


mm 


fl-l  n2 

Unsh  1 
tri 


-1  I  V-''  tx — 1  1 

1  +  Em=l  a 

y1 


J  —  Ti,  .  .  .  ,  T#_i 

j  =  Tr 


and  Ir-\  =  1. 
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Proof.  Denote 


v(y)  = 


(c;y)s 


1  +  llyll2 

Note  that  ||y^||  =  r1.  Since  ri  e  T^y®  is  also  a  feasible  point  of  Subproblem  1.  Thus,  rj  (y^)  < 
77  (y(1))  due  the  optimality  of  yn^  in  Subproblem  1.  This  means  that  there  is  no  need  to  consider 
Subproblem  0.  Therefore,  to  solve  the  optimization  problem,  we  only  need  to  check y (1) , . . . ,  y!/;  ~lj 


and  find  the  one  with  the  largest  77  (y  W) .  However,  for  i  —  1, . . . ,  R  —  2,  if 

V  7  l^m= 1  ° 

which  is  equivalent  to  f  =  0, 


1  +  Em  =  l  arm  >  1-1 


Oj  j  T\ ,  .  .  .  ,  Tj 

Cj  j  =  Ti+i,...,TR. 


and  _ 

i 

^rf+1  ||Cri+i,...,Ti{  ||  +  ^  ^  \aTj\"  =  ri+ 1- 

3= 1 

Since  ri+1  e  Ti+1,  yW  is  a  feasible  point  of  Subproblem  i  +  1.  Thus,  77  (yW)  <  77  (y^+1^)  due  the 
optimality  of  ya+1)  in  Subproblem  7  +  1.  This  means  that  there  is  no  need  to  consider  Subproblem 
i  +  l.  Thus,  we  we  only  need  to  check  those  y^’s  with  f  =  1.  From  (12),  (15),  (16),  and 
/,  =  1,  we  have  x(?l  =  A~T yW.  To  maximize  the  receive  SNR,  we  only  need  to  check  those  x^’s 
for  7  =  1, . . . ,  R  —  1  with  I7  —  1  and  find  the  one  with  the  largest  receive  SNR.  Thus,  (14)  is 
obtained.  □ 


3.2  Discussion 

It  is  natural  to  expect  the  power  control  at  relays  to  undergoe  an  on-or-off  scenario:  a  relay  uses  its 
maximum  power  if  its  channels  are  good  enough  and  otherwise  not  to  cooperate  at  all.  Our  result 
shows  otherwise.  The  optimal  power  used  at  a  relay  can  be  any  value  between  0  and  its  maximal 
power.  In  many  situations,  a  relay  should  use  partial  of  its  power,  whose  value  is  determined  not 
only  by  its  own  channels  but  all  others’  as  well.  This  is  because  every  relay  has  two  effects  on  the 
transmission.  For  one,  it  helps  the  transmission  by  forwarding  the  information,  while  for  the  other, 
it  harms  the  transmission  by  forwarding  noise  as  well.  Its  transmit  power  has  a  non-linear  effect  on 
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the  powers  of  both  the  signal  and  the  noise,  which  makes  the  optimization  solution  not  an  on-or-off 
one,  not  a  decoupled  one,  and,  in  general,  not  even  a  differentiable  function  of  channel  coefficients. 


As  shown  in  Theorem  1  and  Lemma  3,  if  x*  =  for  some  i  E  [1,  R  —  1],  the  fraction  of  power 


used  at  relay  j  satisfies  Gy  =  1  for  j  =  T\, ... .  n  and  ay  =  dp)j  for  j  =  t1+  \  . . . . ,  tr,  where 


for  i  —  1, . . . ,  R  —  2  and 


Thus,  the  i  relays  whose  0’s  are  the  largest  use  their  maximal  powers.  Since  i  >  1,  there  is  at  least 
one  relay  that  uses  its  maximum  power.  This  tells  us  that  the  relay  with  the  largest  0  always  uses  its 
maximal  power.  The  remaining  R  —  i  relays  whose  0’s  are  smaller  only  use  parts  of  their  powers. 
For  j  =  Ti+ 1, . . . ,  Tr,  the  power  used  at  the  jth  relay  is  ajPj  =  dj^Pj  =  d-\fj/gj\2  (1  +  \fj\2P0), 
which  is  proportional  to  | fj/gj\2  (1  +  \fj\2P0)  since  di  is  a  constant  for  each  channel  realization. 
Although  Pj  does  not  appear  explicitly  in  the  formula,  it  affects  the  decision  of  whether  the  jth  relay 
should  use  its  maximal  power.  Actually,  in  determining  whether  a  relay  should  use  its  maximal 
power,  not  only  do  the  channel  coefficients  and  power  constraint  at  this  relay  account,  but  also  all 
other  channel  coefficients  and  power  constraints.  The  power  constraint  of  the  transmitter,  PQ,  plays 
a  roll  as  well. 

Due  to  these  special  properties  of  the  optimal  power  control  solution,  it  can  be  implemented 
distributively  with  each  relay  knowning  only  its  own  channel  information.  In  the  following,  we 
propose  two  distributed  strategies.  One  is  for  networks  with  a  small  number  of  relays,  and  the  other 
is  more  economical  in  networks  with  a  large  number  of  relays. 

The  receiver,  which  knows  all  channels,  can  solve  the  power  control  problem.  When  the  number 
of  relays,  R,  is  small,  the  receiver  broadcases  the  indexes  of  the  relays  that  use  their  full  powers  and 
the  coefficient  di.  If  relay  j  hears  its  own  index  from  the  receiver,  it  will  use  its  maximal  power 
to  transmit  during  the  second  step.  Otherwise,  it  will  use  power  d2\fj/gj\2  (1  +  \fj\2P0).  The  bits 
needed  for  the  feedback  is 


i  log  R  +  Bi  <  R  log  R  +  Bi , 
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where  i  is  the  number  of  relays  that  use  their  maximal  powers  and  B\  is  the  number  of  bits  needed 
in  broadcasting  the  real  number  d*.  Instead,  the  receiver  can  also  broadcast  two  real  numbers:  d,  and 
a  real  number  d  that  satisfies  (j)Ti  >  d  >  (f>Ti+1 .  Relay  j  calculates  its  own  (f)j.  If  (f>j  >  d,  relay  j  uses 
its  maximal  power.  Otherwise,  it  uses  power  d2\fj/gj\ 2  (1  +  \fj\2Po).  The  number  of  bits  needed 
for  the  feedback  is  2 B]_.  Thus,  when  R  is  large,  this  strategy  needs  less  bits  of  feedback  compared 
to  the  first  one. 

Networks  with  an  aggregate  power  constraint  P  on  relays  were  analyzed  in  [24].  In  this  case, 
with  the  same  notation  in  Section  3.1,  Pj  =  P  and  i  aj  —  1-  The  optimal  solution  is 

I  fj9*\y/i+\fipPb 

_  \fi\2Po+\9j\2P+l 

CX,  j  /  • 

/vT R  |/mgm|2(l  +  |/m|2Po) 

V  ^rn=l  (|/m|2Po  +  |gm|2P+l)2 

Qj  is  a  function  of  its  own  channels  fj ,  gj  only  and  an  extra  coefficient  c  =  \jYlm= i  (I’^J^+l^'i^p+i)2  ’ 
which  is  the  same  for  all  relays.  Therefore,  this  power  allocation  can  be  done  distributively  with  the 
extra  knowledge  of  one  single  coefficient  c,  which  can  be  broadcasted  by  the  receiver.  In  our  case, 
every  relay  has  a  separate  power  constraint.  This  is  a  more  practical  assumption  in  sensor  networks 
since  every  sensor  or  wireless  device  has  its  own  battery  power  limit.  The  power  control  solutions 
of  the  two  cases  are  totally  different. 

If  relay  selection  is  used  and  only  one  relay  is  allowed  to  cooperate,  it  can  be  proved  easily  that 
we  should  choose  the  relay  with  the  highest 


i)  fail 2 

i  +  i/^o  +  y2iy 


We  call  h  the  relay  selection  function  since  a  relay  with  a  larger  hj  results  in  a  higher  receive  SNR. 
While  all  relays  are  allowed  to  cooperate,  the  concepts  of  the  best  relay  and  relay  selection  function 
are  not  clear.  Since  the  power  control  problem  is  a  coupled  one,  it  is  hard  to  measure  how  much 
contribution  a  relay  has.  As  discussed  before,  in  network  beamforming,  a  relay  with  a  larger  o3  does 
not  necessarily  use  a  larger  power  or  has  more  contribution.  But  we  can  conclude  that  if  <f>k  >  fa, 
the  fraction  of  power  used  at  relay  k,  ak,  is  no  less  than  the  fraction  of  power  used  at  relay  /,  ai.  It 
is  worth  to  mention  that  in  network  beamforming,  relays  with  larger  enough  fas  use  their  maximal 
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powers  no  matter  what  their  maximal  powers  are.  Actually,  it  is  not  hard  to  see  that  if  at  one  time 
channels  of  all  relays  are  good,  every  relay  should  use  its  maximum  power. 

3.3  Simulation  Results 

In  this  section,  we  show  simulated  performance  of  network  beamforming  and  compare  it  with  per¬ 
formance  of  other  existing  schemes.  Figures  2(a)  and  2(b)  show  performance  of  networks  with 
Rayleigh  fading  channels  and  the  same  power  constraint  on  the  transmitter  and  relays.  In  other 
words,  fi,gi  are  CAf(0, 1)  and  P0  =  P1  =  ■  ■  ■  =  PR  =  P.  The  horizontal  axis  of  the  figures 
indicates  P.  In  Fig.  2(a),  simulated  block  error  rates  of  network  beamforming  with  optimal  power 
control  are  compared  to  those  of  best-relay  selection,  Larsson’s  scheme  in  [24]  with  total  relay 
power  P,  distributed  space-time  coding  in  [10],  and  amplify-and-forward  without  power  control 
(every  relay  uses  its  maximal  power)  in  a  2-relay  network.  The  information  symbol  s  is  modulated 
as  BPSK.  We  can  see  that  network  beamforming  with  optimal  power  control  outperforms  all  other 
schemes.  It  is  about  0.5dB  and  5dB  better  than  Larsson’s  scheme  and  best-relay  selection,  respec¬ 
tively.  With  perfect  channel  knowledge  at  relays,  it  is  7dB  better  than  Alamouti  distributed  space- 
time  coding,  which  needs  no  channel  information  at  relays.  Amplify-and-forward  with  no  power 
control  only  achieves  diversity  1,  best-relay  selection  and  distributed  space-time  coding  achieve  a 
diversity  slightly  less  than  two,  while  network  beamforming  and  Larsson’s  scheme  achieve  diversity 
2.  Fig.  2(b)  shows  simulated  performance  of  a  3-relay  network  under  different  schemes.  Similar 
diversity  results  are  obtained.  But  for  the  3-relay  case,  network  beamforming  is  about  1.5dB  and 
6dB  better  than  Larsson’s  scheme  and  best-relay  selection,  respectively. 

In  Fig.  3(a),  we  show  performance  of  a  2-relay  network  in  which  P0  =  P\  =  P  and  P2  =  P/2. 
That  is,  the  transmitter  and  the  first  relay  have  the  same  power  constraint  while  the  second  relay  has 
only  half  the  power  of  the  first  relay.  The  channels  are  assumed  to  be  Rayleigh  fading  channels.  In 
Fig.  3(b),  we  show  performance  of  a  2-relay  network  whose  channels  have  both  fading  and  path-loss 
effects.  We  assume  that  the  distance  between  the  first  relay  and  the  transmitter/receiver  is  1,  while 
the  distance  between  the  second  relay  and  the  transmitter/receiver  is  2.  The  path-loss  exponent  [38] 
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is  assume  to  be  2.  We  also  assume  that  the  transmitter  and  relays  have  the  same  power  constraint, 
i.e.,  P0  =  Pi  =  P2  =  P.  In  both  cases,  distributed  space-time  coding  does  not  apply,  and  Larsson’s 
scheme  applies  for  the  second  case  only.  So,  we  compare  network  beamforming  with  best-relay 
selection  and  amplify-and-forward  with  no  power  control  only.  Performance  of  Larsson’s  scheme 
is  shown  in  Fig.  3(b)  as  well.  Both  figures  show  the  superiority  of  network  beamfoming  to  other 
schemes. 


4  Networks  with  a  Direct  Link 

The  previous  section  is  on  power  control  of  relay  networks  with  no  DL  between  the  transmitter  and 
receiver.  In  this  section,  we  discuss  networks  with  a  DL.  As  in  [8],  there  are  several  scenarios,  which 
we  discuss  separately. 


4.1  Direct  Link  During  the  First  Step  Only 


In  this  subsection,  we  consider  relay  networks  with  a  DL  during  the  first  step  only.  This  happens 
when  the  receiver  knows  that  the  transmitter  is  in  vicinity  and  listens  during  the  first  step,  while  the 
transmitter  is  not  aware  of  the  DL  or  is  unwilling  to  do  the  optimization  because  of  its  power  and 
delay  constraints.  It  can  also  happen  when  the  transmitter  is  in  the  listening  or  sleeping  mode  during 
the  second  step. 

In  this  case,  /30  =  0.  From  (1)  and  (3),  the  system  equations  can  be  written  as 


Xl 

x2 


aoy/Pofo 

a°VF°  ^  v/l+^PPo. 


S  + 


Wi 

QilgilVA  e  —j  arg  fjn 

2  +  Li=i  +^i/lPP0e 


Using  maximum  ratio  combining,  the  ML  decoding  is 


arg  max 


R 


X\  —  a0V^ofoS  +  |  1  + 


<Xj\9i\  Pi 

tl1  +  a2o\M2po 


-1 


X2  -  CKoV  Pq  2_^  =s 


Sf  V71  +  «ol/il2-Po 


The  optimization  problem  is  thus  the  maximization  of  the  total  receive  SNR  of  both  transmission 
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branches,  which  equals 


2 


First,  both  terms  in  the  SNR  formula  increas  as  a0  increases.  Thus,  a*0  =  1,  i.e.,  the  transmitter 
should  use  its  maximum  power.  The  SNR  optimization  problem  becomes  the  one  in  Section  3.1,  in 
which  there  is  no  DL.  Therefore,  the  power  control  of  networks  with  a  DL  during  the  first  step  only 
is  exactly  the  same  as  that  of  networks  without  a  DL.  This  result  is  intuitive.  Since  with  a  DL  during 
the  first  step  only,  operations  at  both  the  transmitter  and  relays  keep  the  same  as  networks  without 
the  DL.  The  only  difference  is  that  the  receiver  obtains  some  extra  information  from  the  transmitter 
during  the  first  step,  and  it  can  use  the  information  to  improve  the  performance  without  any  extra 
cost.  For  the  single -relay  case,  it  can  be  proved  easily  that  to  maximize  the  receive  SNR,  the  relay 
should  use  its  maximal  power  as  well,  that  is,  a \  =  1. 

4.2  Direct  Link  During  the  Second  Step  Only 

In  this  subsection,  we  consider  relay  networks  with  a  DL  during  the  second  step  only.  This  happens 
when  the  transmitter  knows  that  the  receiver  is  at  vicinity  and  determines  to  do  more  optimization 
to  allocate  its  power  between  the  two  transmission  steps.  However,  the  receiver  is  unaware  of  the 
DL  and  is  not  listening  during  the  first  step.  It  can  also  happen  when  the  receiver  is  in  transmitting 
or  sleeping  mode  during  the  first  step. 

In  this  case,  X\  =  0  and  x2  is  given  in  (3).  The  receive  SNR  can  be  calculated  to  be 


i  ,  \x\'2Pi 

^  2^i=  1  l+C*g|/j|2P0 

First,  we  show  that  aft  +  /3q  should  take  its  maximal  value  1,  i.e.,  the  transmitter  should  use  all  its 


power.  Assume  that  +  $ g  <  1  is  the  optimal  solution.  Define  T0  =  a/1  —  d^.  Wc  have  T0  >  T0. 
Therefore,  SNR(a0 ,  d0)  <  SNR(a0,  /30).  This  contradicts  the  assumption  that  (d0,  /%)  is  optimal. 
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Define 


\9i\VPi  t  «o \fi9i\y/Pi  -  bi  i  ,.  rM  ,  *  ._i 

Ui  =  21  f  1 2  P  1  ’  bi  =  21  f  12  o  ’  Q  =  F>  ^  =  dlag{a)>  aIld  J/<  =  “i- 

v  1  +  aol/*r-^o  vl  +  aol/*r-Po  a* 


The  receiver  SNR  can  be  calculated  to  be 

(V1  -  «ol/ol  +  (b,x) 
^(a0,x)  =  P0  — 


=  Pn 


V1  -  «ol/ol  +  (c,y> 


1  +  Px||2  ”  i  +  lly|l2 

For  any  fixed  a0,  we  can  optimize  «i, . . . ,  aR  following  the  analysis  in  Section  3.1.  The  following 
theorem  can  be  proved. 

Theorem  2.  For  any  fixed  a0  G  (0, 1),  order  (f  =  'f-  as 


^Tl  ^  ^  ff R- 


The  receive  SNR  is  maximized  at 


x  =  arg  max 


V1  “  «ol/ol  +  (b,xW) 


t(0 


1  +  ||ylx«||2 

where  for  i  —  0, . . . ,  R  —  2,  the  R-dimensional  vector  is  defined  as 


a(»)  _ 
xj  ~ 


i+E 


m= 1  r m 


\A-^ol/ol+Em=l  N 


J  =  Ti,.  .  .,Ti 
J  Ts+lj  •  •  •  j  Tr 


and  the  binary  functions  /,  is  defined  as 


JR-l) 


is  defined  as 


h  = 


0  // 


i+E 


% 


> 


i-l 


\A^!l/o|+EUiNm  “  fi+1 

1  otherwise 


x 


(R-l) 


mm 


i+E 


R-l 

m=l  afm 


\/1_"ol/o|  +  Em=l 


yTji 


,1 


and  /R_!  =  1. 


J  h)  ■  ■  ■  i  D?—  i 

3  =  Tr 
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Proof.  The  proof  of  this  theorem  follows  the  one  of  Theorem  1  and  the  lemmas  it  uses. 


□ 


As  discussed  in  Section  3.1,  for  networks  with  no  DL,  there  is  no  need  to  consider  the  solution 
of  Subproblem  0.  Here  it  is  different.  Define  f\  =  071||c||.  If  we  denote  the  solution  of  Subproblem 
0,  max^igp^o^y^a  (a1tf^2>  »  as  y(0),  because  of  the  existance  of  the  DL  during  the  second  step, 
it  is  possible  that  ||y(°)  ||  <  f  Thus,  we  need  to  check  one  extra  vector  x^0)  =  A~ 7 y(°)  and  binary 
function  I0. 

Now  we  discuss  the  optimization  of  a0.  We  first  consider  the  case  of  a0  6  (0, 1).  For  any  given 


x  = 


oiy  ■  ■  olr 


,  the  O'q  that  maximizes  the  receive  SNR  satisfies  Jyy  =  0.  Thus,  the  optimal 


can  be  found  numerically  by  solving  =  0.  It  can  be  proved  easily  that  >  0  when  a0  — >  0+ 


dip 


8a0 


da0 


and  <  0  when  cto  1  .  Thus,  the  maximum  of  f  is  reached  inside  (0, 1). 


When  the  power  at  the  transmitter  is  high  (P0  1),  the  receive  SNR  can  be  approximated  by 

Pq  (y/T^M  + 


■f(a0,x)  «  d(a0)  = 


1  +  d2/a.Q 


where  d i  =  YfLi  ai\di\  '/P  and  d2  =  ^  Jf  ,=\  at\fh/ fi\2Pf  It  can  be  calculated  straightfor¬ 
wardly  that  for  a0  G  (0, 1), 


4Po  ^\/l  —  Qtol/ol  +  d\  j 

da0  a^y/T^of(l  +  d2/al)2  _ 


dd 


1/oK  —  26|/0|aQ  +  b\f0\  +  ab\J~\ 


a? 


and 


dd 

dan 


=  0  |/o|  a0  -  M2\fo\  a0  +  2d2\f0\  a0  +  d2(4|/0|  -  df)aQ  +  d2{d1  -  |/0|  )  =  0 


This  is  a  quartic  equation  of  olq,  whose  solutions  can  be  calculated  analytically.  Note  that  yy  >  0 
when  a0  — >  0+  and  <  o  when  a0  — >•  1_.  Thus  the  maximum  of  d  is  reached  inside  (0, 1).  An 
approximate  solution  of  a0  can  thus  be  obtained  analytically  at  high  transmit  powers. 

Now  we  consider  the  cases  of  a0  =  0  and  a0  =  1.  If  a0  =  0,  the  system  degrades  to  a  point- 
to-point  one  since  only  the  DL  works.  Thus,  the  receive  SNR  is  |/0|2Po.  For  o0  =  1,  we  can 
obtain  the  optimal  x  using  Theorem  2.  Thus,  we  obtain  three  sets  of  a  and  x  for  the  three  cases: 
«o  G  (0, 1),  a0  =  0,  and  a0  =  1,  respectively.  The  optimal  solution  of  the  system  is  the  set  of  a0 
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and  x  corresponding  to  the  largest  receive  SNR.  The  power  control  problem  in  networks  with  a  DL 
during  the  second  step  only  can  thus  be  solved  using  the  following  recursive  algorithm. 


Algorithm  1. 


1.  Initialization:  Set  yR,remous'>  —  \R>  the  R-dimensional  vector  of  all  ones,  S N R(J>remous'>  —  q, 
and  count  =  0.  Set  the  maximal  number  of  iterations  iter  and  the  threshold  thre. 

2.  Optimize  oq  with  x  =  x'/'" 1  lous\  Denote  the  solution  as  Oq.  We  can  either  do  this  numeri¬ 
cally  or  calculate  the  high  SNR  approximation. 


3.  With  a0  =  a^,  find  the  x  that  maximizes  the  receive  SNR  using  Theorem  2.  Denote  it  as  X] . 
Calculate  SNRX  =  ^(aj^xx). 


4.  Set  count  =  count+l.  If  count  <  iter  and 


SNRi  -  SNR 


( previous ) 
1 


>  thre,  set  x^ 


(previous) 


xi,  SN R(frf  Vlous>  =  SNRi,  and  go  to  step  2. 


5.  Find  the  solution  o/x  with  o0  =  1  using  Theorem  2.  Denote  this  solution  as  x2. 

6.  The  optimal  solution  is:  («q,  x*)  =  argmax  <  ^(a^,  Xx),  ^(1,  x2),  ^(0,  0^)  i. 


Similarly,  the  distributive  strategies  proposed  in  Section  3.2  can  be  applied  here. 


4.3  Direct  Link  During  Both  Steps 

In  this  subsection,  we  consider  relay  networks  with  a  DL  during  both  the  first  and  the  second  steps. 
This  happens  when  both  the  transmitter  and  the  receiver  know  that  they  are  not  too  far  away  from 
each  other  and  decide  to  communicate  during  both  steps  with  the  help  of  relays  during  the  second 
step. 

From  (1)  and  (3),  the  system  equation  can  be  written  as 

Xi 

X2 


C^os/Pofo 


w2  +  z2i= i 


W\ 

aAgiWPj  c  — 7  arg.fi r 


vT+agl  fi\2Po 
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Similar  to  the  networks  discussed  in  Section  4.1,  the  maximum  ratio  combining  results  in  the 
following  ML  decoding: 


arg  max 

S 


X\  —  O'  o 


y/Poft 


o  s 


+ 


^  ai\9i\2Pi 


-1 


l+^WPo, 

The  total  receive  SNR  of  both  transmission  branches  can  be  calculated  to  be 

2 


«o^o|/o|2  + 


1  ai\9i\2pi 

"l"  ^i=1  i+«gl/il2Ro 


=  a20P0\fo\2  +  ip(ao,x). 


The  same  as  the  networks  in  Section  4.2,  a g  +  should  take  its  maximal  value,  which  is  1.  That  is, 
A)  =  \/i —  ao-  Similar  to  the  SNR  optimization  in  Section  4.2,  for  any  given  a0  e  (0, 1),  the  SNR 
maximization  is  the  same  as  the  maximization  of  ip,  which  is  solved  by  Theorem  2.  But  due  to  the 
difference  in  the  receive  SNR  formula,  the  optimal  q0  given  an, . . . ,  an  is  different.  It  is  the  solution 
of  2o0/ o|/()|2  +  =  0.  When  the  DL  exists  during  both  steps,  the  case  of  a0  =  0,  whose  receive 

SNR  is  \f0\2Po  will  never  outperform  the  case  of  a0  =  1,  whose  receive  SNR  is  |/0|2P0  +  V’(l) x)  f°r 
some  x.  Thus,  the  case  ao  =  0  needs  not  to  be  considered.  The  power  control  problem  in  networks 
with  a  DL  during  both  steps  can  thus  be  solved  using  the  following  recursive  algorithm. 


Algorithm  2. 

1.  Initialization:  Set  x(Premous')  —  \R>  sNR^>remous'1  =  0,  and  count  =  0.  Set  the  maximal 
number  of  iterations  iter  and  the  threshold  thre. 

2.  Optimize  ao  with  x  =  ■x^rremous'> _  [)enoie  the  solution  as  Oq  \  We  can  do  this  numerically. 

3.  With  ao  =  a:^1 ,  find  the  x  that  maximizes  ip  using  Theorem  2.  Denote  it  as  xt.  Calculate 


SNRi  =  (a^y\fo\2P0  +  ip(^\^). 


4.  Set  count  =  count+1.  If  count  <  iter  and 


SNR x  -  SNR{rmous)  >  thre,  set  = 


xi,  SNR[pit  “OUh)  =  SNR\  and  go  to  step  2. 


5.  Find  the  solution  o/x  with  a0  =  1  using  Theorem  2.  Denote  this  solution  as  x2. 
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6.  The  optimal  solution  is:  (a^,  x*)  =  arg  max  |  \fo\2Po  +  V;(Qo'\  xi);  \fo\2Po  +  ^(1,  x2) 

Again,  the  distributive  strategies  proposed  in  Section  3.2  can  be  applied  here. 

4.4  Performance  Comparison 

In  this  subsection,  we  compare  single -relay  networks  in  which  the  power  constraints  at  the  transmit¬ 
ter  and  the  relay  are  same,  i.e.,  P0  —  Pi  —  P.  The  channels  are  assumed  to  have  both  the  fading 
and  path-loss  effect.  There  are  four  cases:  no  DL,  a  DL  during  the  first  step  only,  a  DL  during  the 
second  step  only,  and  a  DL  during  both  steps. 

In  Fig.  5,  we  compare  networks  in  which  the  distance  of  every  link  is  the  same,  i.e.,  the  three 
nodes  are  vertexes  of  an  equilateral  triangle  with  unit-length  edges  as  shown  in  Fig.  4(a).  We  can  see 
that  the  network  with  no  DL  has  diversity  1  while  networks  with  a  DL  and  power  control  achieve 
diversity  2.  The  network  with  a  DL  during  the  first  step  performs  less  than  0.5dB  better  than  the 
network  with  a  DL  during  the  second  step  only,  while  the  network  with  a  DL  during  both  steps 
performs  the  best  (about  ldB  better  than  the  network  with  a  DL  during  the  first  step  only).  To 
illuminate  the  effect  of  power  control,  we  show  performance  of  networks  whose  transmit  power  at 
the  relay  and  transmitter  are  fixed.  For  the  network  with  a  DL  during  the  first  step  only,  there  is  no 
power  control  problem  since  it  is  optimal  for  both  the  transmitter  and  the  relay  to  use  their  maximal 
powers.  For  the  other  two  cases,  we  let  the  transmitter  uses  half  of  its  power,  P/2,  to  each  of  the  two 
steps  and  the  relay  always  uses  its  maximum  power  P.  We  can  see  that,  if  the  DL  only  exists  during 
the  second  step,  without  power  control,  the  achievable  diversity  is  1.  At  block  error  rates  of  10“  2 
and  10~“3,  it  performs  3  and  6dB  worse,  respectively.  For  networks  with  a  DL  during  both  steps, 
power  control  results  in  a  1.5dB  improvement. 

In  Fig.  6(a)  and  6(b),  we  show  performance  of  line  networks  with  path-loss  exponents  2  and  3 
respectively.  As  shown  in  Fig.  4(b),  the  three  nodes  are  on  a  line  and  the  relay  is  in  the  middle  of 
the  transmitter  and  receiver.  The  distance  between  the  transmitter  and  receiver  is  assumed  to  be  2. 
The  same  phenomenon  as  in  the  equilateral  triangle  networks  can  be  observed.  The  network  with 
a  DL  during  both  steps  performs  the  best  (about  ldB  better  than  the  network  with  a  DL  during  the 
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first  step  only).  The  network  with  a  DL  at  first  step  only  performs  slightly  better  than  the  one  with  a 
DL  during  the  second  step  only.  But  the  difference  is  smaller  than  that  in  Fig.  5.  The  performance 
difference  between  line  networks  with  and  without  DLs  is  smaller  than  those  in  equilateral  triangle 
networks,  and  it  gets  even  smaller  for  larger  path-loss  exponents.  This  is  because  as  the  distance 
between  the  transmitter  and  receiver  or  the  path-loss  exponent  is  larger,  the  quality  of  the  DL  is 
lower.  Therefore,  the  improvement  due  to  this  link  is  smaller.  For  both  cases,  power  control  results 
in  a  1.5dB  improvement  when  the  DL  link  exists  for  both  steps  and  a  higher  diversity  when  the  DL 
exists  for  the  second  step  only. 

Then  we  work  on  the  random  network  in  Fig.  4(c),  in  which  the  relay  locates  randomly  and 
uniformly  within  a  circle  in  the  middle  of  the  transmitter  and  the  receiver.  The  distance  between  the 
transmitter  and  the  receiver  is  assumed  to  be  2.  The  radius  of  the  circle  is  denoted  as  r.  We  assume 
that  r  <  1.  This  is  a  reasonable  model  for  ad  hoc  wireless  networks  since  if  communications  between 
two  nodes  is  allowed  to  be  helped  by  one  other  relay,  one  should  choose  a  relay  that  is  around  the 
middle  of  the  two  nodes.  In  other  words,  the  distance  between  the  relay  and  the  transmitter  or 
receiver  should  be  shorter  than  that  between  the  transmitter  and  receiver. 

We  work  out  the  geometry  first.  As  in  Fig.  4(c),  we  denote  the  positions  of  the  transmitter, 
the  receiver,  the  relay,  and  the  middle  point  of  the  transmitter  and  the  receiver  as  A,C,  D,  and  B, 
respectively.  Denote  the  angle  of  AB  and  BD  as  6  and  the  length  of  BD  as  p.  The  lengths  of  AD 
and  CD  are  thus  \J  1  +  p2  —  2pcos0  and  \Jl  +  p2  +  2pcos  6.  Since  D  is  uniformly  distributed 
within  the  circle,  6  is  uniform  in  [0, 7r)  and  the  pdf  and  cdf  of  p  can  be  calculated  to  be 


Ip  x 

p(p)  =  —  and  P(p<x)  =  — , 


respectively.  Define  Y  =  r\[X.  If  X  is  uniform  on  (0, 1),  it  can  be  proved  that 


Thus,  Y  has  the  same  distribution  as  p.  Therefore,  we  generate  Y  to  represent  p. 

Fig.  7  shows  performance  of  random  networks  with  path-loss  exponent  2  and  r  =  1/2.  We  can 
see  that  the  same  phenomenon  as  in  line  networks  can  be  observed.  With  a  DL  at  both  steps,  the 
random  network  performs  about  ldB  worse  than  the  line  network. 
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5  Conclusions  and  Future  Work 


In  this  paper,  we  propose  the  novel  idea  of  beamforming  in  wireless  relay  networks  to  achieve 
both  diversity  and  array  gain.  The  scheme  is  based  on  a  two-step  amplify-and-forward  protocol. 
We  assume  that  each  relay  knows  its  own  channels  perfectly.  Unlike  previous  works  in  network 
diversity,  the  scheme  developed  here  uses  not  only  the  channels’  phase  information  but  also  their 
magnitude.  Match  filters  are  applied  at  the  transmitter  and  relays  during  the  second  step  to  cancel  the 
channel  phase  effect  and  thus  form  a  coherent  beam  at  the  receiver,  in  the  mean  while,  optimal  power 
control  is  performed  based  on  the  channel  magnitude  to  decide  the  power  used  at  the  transmitter  and 
relays.  The  power  control  problem  for  networks  with  any  numbers  of  relays  and  no  direct  link  is 
solved  analytically.  The  solution  can  be  obtained  with  a  complexity  that  is  linear  in  the  number  of 
relays.  The  power  used  at  a  relay  depends  on  not  only  its  own  channels  nonlinearly  but  also  all  other 
channels  in  the  network.  In  general,  it  is  not  even  a  differentiable  function  of  channel  coefficients. 
Simulation  with  Rayleigh  fading  and  path-loss  channels  show  that  network  beamforming  achieves 
the  maximum  diversity  while  amplify-and-forward  without  power  control  achieves  diversity  1  only. 
Network  beamforming  also  outperformes  other  cooperative  strategies.  For  example,  it  is  about  4dB 
better  than  best-relay  selection. 

Relay  networks  with  a  direct  link  between  the  transmitter  and  receiver  are  also  considered  in  this 
paper.  For  networks  with  a  direct  link  during  the  first  step  only,  the  power  control  at  relays  and  the 
transmitter  is  exactly  the  same  as  that  of  networks  with  no  direct  link.  For  networks  with  a  direct 
link  during  the  second  step  only  and  networks  with  a  direct  link  during  both  steps,  the  solutions  are 
different.  Recursive  numerical  algorithms  for  the  power  control  at  both  the  transmitter  and  relays  are 
given.  Simulated  performance  of  single-relay  networks  with  different  topologies  shows  that  optimal 
power  control  results  in  about  1.5dB  improvement  in  networks  with  a  direct  link  at  both  steps  and  a 
higher  diversity  in  networks  with  a  direct  link  at  the  second  step  only. 

We  have  just  scratched  the  surface  of  a  brand-new  area.  There  are  a  lot  of  ways  to  extend  and 
generalize  this  work.  First,  it  is  assumed  in  this  work  that  relays  and  sometimes  the  transmitter 
know  their  channels  perfectly,  which  is  not  practical  in  many  networks.  Network  beamforming 
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with  limited  and  delayed  feedback  from  the  receiver  is  an  important  issue.  In  multiple- antenna 
systems,  beamforming  with  limited  and  delayed  channel  information  feedback  has  been  widely 
probed.  However,  beamforming  in  networks  differs  from  beamforming  in  multiple-antenna  systems 
in  a  couple  of  ways.  In  networks,  it  is  difficult  for  relays  to  cooperate  while  in  a  multiple- antenna 
system,  different  antennas  of  the  transmitter  can  cooperate  fully.  There  are  two  transmission  steps  in 
relay  networks  while  only  one  in  multiple-antenna  systems,  which  leads  to  different  error  rate  and 
capacity  calculation  and  thus  different  designs.  Second,  the  relay  network  probed  in  this  paper  has 
only  one  pair  of  transmitter  and  receiver.  When  there  are  multiple  transmitter-and-receiver  pairs,  an 
interesting  problem  is  how  relays  should  allocate  their  powers  to  aid  different  communication  tasks. 
Finally,  the  two-hop  protocol  can  be  generalized  as  well.  For  a  given  network  topology,  one  relevant 
question  is  how  many  hops  should  be  taken  to  optimize  the  criterion  at  consideration,  for  example, 
error  rate  or  capacity. 
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Figure  1:  Wireless  relay  network. 


(a)  2-relay  network  (b)  3-relay  network 

Figure  2:  Networks  with  fading  channels  and  same  power  constraint  for  all  nodes. 
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(a)  Network  with  different  relay  powers  (b)  Network  with  path-loss  plus  fading  channels 

Figure  3:  2-relay  networks  with  different  relay  powers  and  pass-loss  plus  fading  channels. 


Relay  Region 


(a)  Triangle  network 


(b)  Line  network 


(c)  Random  network 


Figure  4:  Network  topology. 
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Figure  5:  Equilateral  triangle  network. 


(a)  Path-loss  exponent  2  (b)  Path-loss  exponent  3 


Figure  6:  Single-relay  line  network. 
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Figure  7:  Single-relay  random  network  with  pass-loss  exponent  2  and  r  =  1/2. 
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